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Abstract 

Electron transport properties through single conjugated molecules sandwiched between two non- 
superconducting electrodes are studied by the use of Green's function technique. Based on the tight- 
binding model, we do parametric calculations to characterize the electron transport through such molecu- 
lar bridges. The electron transport properties are significantly influenced by (a) the existence of localizing 
groups in these conjugated molecules and (b) the molecule to electrode coupling strength, and, here we 
focus our results in these two aspects. 
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1 Introduction 

Molecular electronics and transport have attracted 
much more attention since molecules constitute 
promising building blocks for future generation of 
nanoelectronic devices. Following experimental de- 
velopments, theory can play a major role in un- 
derstanding the new mechanisms of conductance. 
The single-molecule electronics plays a significant 
role in designing and developing the future nano- 
electronic circuits, but, the goal of developing a re- 
liable molecular-electronics technology is still over 
the horizon and many key problems, such as device 
stability, reproducibility and the control of single- 
molecule transport need to be solved. Electronic 
transport through molecules was first studied theo- 
retically in 1974 {![. Then lot of experiments 
[51 [S| have been performed through molecules placed 
between two metallic electrodes with few nanome- 
ter separation. The operation of such two-terminal 
devices is due to an applied bias. Current passing 
across the junction is strongly nonlinear function 
of applied bias voltage and its detailed description 
is a very complex problem. The complete knowl- 
edge of the conduction mechanism in this scale is 
not well understood even today. In many molec- 
ular devices, electronic transport is dominated by 
conduction through broadened HOMO or LUMO 
states. In contrast here we find that the transport 
through single conjugated molecules can be con- 
trolled very sensitively by introducing the localizing 
groups in these molecules. This sensitivity opens up 
new possibilities for novel single-molecule sensors. 
Electron conduction through molecules strongly de- 
pends on (a) the delocalization of the molecular 
electronic orbitals and (b) their coupling strength 
to the two electrodes. In a very recent experi- 
ment. Tali Dadosh et al. [2] have measured conduc- 
tance of single conjugated molecules and predicted 
that the existence of localizing groups in a conju- 
gated molecule suppresses the electrical conduction 
through the molecule. These results motivate us 
to study the electron transport through such con- 
jugated molecules. 

The aim of the present paper is to reproduce an 
analytic approach based on the tight-binding model 
to investigate the electron transport properties for 
the model of single conjugated molecules taken in 
their experiment [2]. Several ab initio methods are 
used for the calculation of conductance [3 [51 [SI [TOl 
[TTl [T2] , yet it is needed the simple parametric ap- 
proaches [I3l[l4l[15l[l6l[l2l[l8l[19l[20l[2l]for this 
calculation. The parametric study is motivated by 



the fact that it is much more flexible than that of 
the ab initio theories since the later theories are 
computationally very expensive and here we focus 
our attention on the qualitative effects rather than 
the quantitative ones. This is why we restrict our 
calculations on the simple analytical formulation of 
the transport problem. 

The scheme of the paper is as follow. In Section 2, 
we give a very brief description for the calculation 
of transmission probability and current through a 
finite size conductor sandwiched between two one- 
dimensional (ID) metallic electrodes. Section 3 fo- 
cuses the results of conductance-energy {g-E) and 
current-voltage {I-V) characteristics for the single 
conjugated molecules and study the effects of lo- 
calizing groups in the above mentioned quantities. 
Finally, we summarize our results in Section 4. 

2 A glimpse onto the theoret- 
ical formulation 

Here we describe very briefly about the methodol- 
ogy for the calculation of transmission probability 
(T), conductance {g) and current (/) through a fi- 
nite size conducting system attached to two semi- 
infinite metallic electrodes by using the Green's 
function technique. 

Let us first consider a ID conductor with N num- 
ber of atomic sites (array of filled circles) connected 
to two semi-infinite electrodes, namely, source and 
drain, as presented in Fig. [1] The conducting sys- 
tem in between the two electrodes can be an array 
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Figure 1: Schematic view of a ID conductor with 
N number of atomic sites (filled circles) attached to 
two electrodes through the sites 1 and N, respec- 
tively. 

of few quantum dots, or a single molecule, or an 
array of few molecules, etc. At low voltages and 
temperatures, the conductance of the conductor can 
be written by using the Landauer conductance for- 
mula. 



where g is the conductance and T is the transmis- 
sion probability of an electron through the conduc- 
tor. The transmission probability can be expressed 



2 



in terms of the Green's function of the conductor 
and the coupHng of the conductor to the two elec- 
trodes by the expression, 



T = Tr [TsG^cTdG^c] 



(2) 



where Gq and Gq are respectively the retarded and 
advanced Green's function of the conductor. Ts 
and r £) are the coupling terms of the conductor due 
to the coupling to the source and drain, respectively. 
For the complete system, i.e., the conductor and the 
two electrodes, the Green's function is defined as, 



(3) 



where e = E + irj. E is the injecting energy of the 
source electron and is a very small number which 
can be put as zero in the limiting approximation. 
The above Green's function corresponds to the in- 
version of an infinite matrix which consists of the 
finite conductor and two semi-infinite electrodes. It 
can be partitioned into different sub-matrices those 
correspond to the individual sub-systems. 

The effective Green's function for the conductor 
can be written as, 



Gc 



(4) 



where He is the Hamiltonian for the conductor 
sandwiched between the two electrodes. The single 
band tight-binding Hamiltonian for the conductor 
within the non-interacting picture can be written in 
the following form. 
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where c| (q) is the creation (annihilation) opera- 
tor of an electron at site i, e^'s are the site ener- 
gies and t is the nearest-neighbor hopping integral. 
Here Eg = hg^gshsc and Ed = hDcgoh^DC ^"^'^ 
the self-energy terms due to the two electrodes, gs 
and go are respectively the Green's function for the 
source and drain, hsc and hue are the coupling 
matrices and they will be non-zero only for the ad- 
jacent points in the conductor, 1 and N as shown 
in Fig. [1] and the electrodes respectively. The cou- 
pling terms Ts and F/j for the conductor can be 
calculated through the expression. 
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where E^^ ^-^ and E°g j-^y are the retarded and ad- 
vanced self-energies, respectively, and they are eon- 
jugate to each other. Datta et al. [H] have shown 
that the self-energies can be expressed like. 
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(7) 



where A|5 £)-j are the real parts of the self-energies 
which correspond to the shift of the energy eigen- 
states of the conductor and the imaginary parts 
^{S,D} of the self-energies represent the broaden- 
ing of these energy levels. This broadening is much 
larger than the thermal broadening and this is why 
we restrict our all calculations only at absolute zero 
temperature. By doing some simple algebra these 
real and imaginary parts of self-energies can also be 
determined in terms of coupfing strength {t[s,d}) 
between the conductor and two electrodes, injec- 
tion energy (E) of the transmitting electron, site 
energy (eg) of the electrodes and hopping strength 
(w) between nearest-neighbor sites in the electrodes. 
Thus the coupling terms Fg and To can be written 
in terms of the retarded self-energy as. 



^{s,D} = -2Im 



^{S,D} 



(8) 



Now all the information regarding the conductor to 
electrode coupling is included into the two self ener- 
gies as stated above and is analyzed through the use 
of Newns- Anderson chemisorption theory [T31 E] . 
The detailed description of this theory is obtained 
in these two references. 

By calculating the self-energies, the coupling 
terms Ts and To can be easily obtained and then 
the transmission probability T can be computed 
from the expression as mentioned in Eq. [21 

Since the coupling matrices hsc and hoc are 
non-zero only for the adjacent points in the conduc- 
tor, 1 and N as shown in Fig. [1] the transmission 
probability becomes, 

TiE, V) = 4Afi(£;, V)A'^j^iE, V)\Gin{E, V)f 

(9) 

The current passing through the conductor is de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current-voltage relation is evaluated from the fol- 
lowing expression [53], 



Ep+eV/2 



m = ^ 
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T{E, V)dE 



(10) 



EF-eV/2 



where Ep is the equilibrium Fermi energy. For the 
sake of simplicity, here we assume that the entire 
voltage is dropped across the conductor-electrode 
interfaces and this assumption does not significantly 
change the qualitative behaviors of the I-V char- 
acteristics. Using the expression of T{E, V) as in 
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Eq. [9]thc final form of I{V) becomes, 

EF+eV/2 
Ep-eV/2 

X \GiN{E,V)\^dE (11) 

Eq.[Tl Eq.[5]and Eq.[Tl]are the final working formu- 
lae for the calculation of conductance g and current- 
voltage characteristics, respectively, for any finite 
size conductor sandwiched between two electrodes. 

With the help of the above formulation, wc shall 
describe the electron transport properties through 
some conjugated molecules (Fig. [5]). For the sake 
of simplicity throughout this article wc use the unit 
c = e = h = 1. 

3 Results and discussion 

This section focuses the conductance-energy (g-E) 
and current-voltage (I-V) characteristics of three 
short single conjugated molecules. These molecules 
are specified as: 1,4-benzenedimethanethiol 



SH 




BPE 



Figure 2: Structures of the three molecules: 
1,4-benzcnedimcthancthiol (BDMT), 4,4'- 
biphenyldithiol (BPD) and bis-(4-mercaptophenyl)- 
ether (BPE) those are attached to two electrodes 
by thiol (S-H) groups. 

(BDMT), in which the molecular conjugation 
is broken near the contacts by a methylene 
group; 4,4'-biphenyldithiol (BPD), a fully conju- 
gated molecule; and bis-(4-mercaptophenyl)-ether 
(BPE), where the molecular conjugation is broken 
by an oxygen atom at the center. The schematic 
representations of these three molecules, with thiol 



groups at the two extreme ends of each molecules, 
are shown in Fig. [2l These molecules are contacted 
to the two semi-infinite ID electrodes by thiol 
(S-H) groups via single channels (same as shown 
schematically in Fig. [T]). In actual experimental ar- 
rangement, two electrodes are constructed by using 
gold (Au) substance and molecule attached to the 
electrodes by thiol (S-H) groups in the chcmisorp- 
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Figure 3: Conductance 5 as a function of the inject- 
ing electron energy E in the weak-coupling limit, 
where (a), (b) and (c) are respectively for the 
BDMT, BPD and BPE molecules. 

tion technique where hydrogen (H) atoms remove 
and sulfur (S) atoms reside. The electron transport 
through such conjugated molecules significantly 
influenced by the presence of localizing groups 
in the molecules and the molecule-to-electrode 
coupling strength. Here, we shall investigate our 
results in these aspects. Throughout the article, 
we discuss the results in two limiting regimes 
depending on the coupling strength of the molecule 
to the electrodes. One is defined as t^s,d} « t, 
the so-called weak-coupling limit. The other one is 
'^{s,D} ~ ^5 the so-called strong-coupling limit. The 
parameters ts and td correspond to the coupling of 
the molecules to the source and drain, respectively 
The values of the different parameters used in 
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our calculations in these two limiting regimes are 
assigned as: T5 = r/j = 0.5; i = 2.5 (weak coupling) 
and Ts = Tu = 2] t = 2.5 (strong-coupling). For 
the side attached electrodes the on-site energy (eg) 
and the nearest-neighbor hopping strength {v) are 
fixed to and 4, respectively. The Fermi energy 
Ep is set at 0. 

In Fig. [3] we display conductance {g) as function 
of injecting electron energy {E) for the three molec- 
ular bridge systems in the limit of weak molecular 
couphng. Figures (HJa), (b) and (c) correspond to 
the results for the bridges with BDMT, BPD and 
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Figure 4: Conductance 17 as a function of the in- 
jecting electron energy E in the strong molecule- 
to-electrode coupling limit, where (a), (b) and (c) 
are respectively for the BDMT, BPD and BPE 
molecules. 

BPE molecules, respectively. Conductance shows 
very sharp resonant peaks for some particular en- 
ergy values, while in almost all other cases it drops 
to zero. These resonant peaks are associated with 
the energy eigenstates of the individual molecules 
that bridges the two reservoirs. Therefore, the con- 
ductance spectrum manifests itself the electronic 
structure of the molecule. At resonances, the con- 
ductance [g) achieves the value 2, and accordingly, 
the transmission probability (T) goes to unity since 



we have the relation g = 2T from the Landauer 
conductance formula with e = /i = 1 in our present 
formulation. For the bridges with BDMT and BPD 
molecules, we see that the resonant peaks have 
very narrow widths, while for the bridge with BPE 
molecule the width of the peaks is almost zero. 
Thus, fine tuning in the energy scale is necessary to 
get the electron conduction through these bridges, 
specially in the case of BPE molecule for the weak- 
coupling limit. The most significant result is that, 
the BPD molecule conducts electron across the zero 
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Figure 5: Current / as a function of the applied 
bias voltage V in the weak molecule-to-electrode 
coupling, where the solid, dotted and dashed lines 
are respectively for the BPD, BDMT and BPE 
molecules. 

energy value, while the other two conduct beyond 
some critical energy values (see Fig. [3]). Therefore, 
we can tune the electron conduction through the 
molecular bridge in a very controllable way. 

In the strong molecular coupling limit, the res- 
onant peaks in the conductance spectra get sub- 
stantial widths as shown in Fig. [D This enhance- 
ment of the resonant widths is due to the broaden- 
ing of the molecular energy eigenstates caused by 
the coupling of the molecules to the side attached 
electrodes in the strong-coupling limit, where the 
contribution comes from the imaginary parts of the 
self-energies [23] . Though for the molecular bridges 
with BDMT and BPD molecules the resonant peaks 
get substantial widths, but, for the bridge with BPE 
molecule, the increment of the widths is very small. 
For this BPE molecule since the increment of the 
width of the resonant peak across the energy E = 2 
is comparatively higher than for the other energy 
values, we observe only one peak across this energy 
value {E = 2, Fig. [4Kc)). Thus, for this molecular 
bridge electron conduction takes place across a par- 
ticular energy value, while in all other energies no 
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electron conduction takes place. This aspect may 
be used to describe switching action in an electronic 
circuit. 

Thus we see that the electron conduction strongly 
depends on the molecule itself and also on the 
strength of the molecular coupling to the side at- 
tached electrodes. The behavior of electron transfer 
through the molecular junction becomes much more 
clearly observed from the current-voltage charac- 
teristics. Current passing through the molecular 
system is computed from the integration procedure 
of the transmission function T . The nature of the 
transmission function is exactly similar to that of 
the conductance spectrum since g = 2T (from the 
Landauer formula) , differ only in magnitude by the 
factor 2. In Fig. [51 we plot the current-voltage char- 
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Figure 6: Current / as a function of the applied bias 
voltage V in the strong moleculc-to-electrode cou- 
pling, where the solid, dotted and dashed curves 
are respectively for the BPD, BDMT and BPE 
molecules. 

acteristics of these three molecular bridges in the 
weak-coupling limit. The solid, dotted and dashed 
curves correspond to the results for the molecular 
bridges with BPD, BDMT and BPE molecules, re- 
spectively. The current shows staircase-like struc- 
ture with fine steps as a function of the applied bias 
voltage. This is due to the discreteness of molecu- 
lar resonances as shown in Fig. [3l With the in- 
crease of the bias voltage, the electrochemical po- 
tentials on the electrodes are shifted and eventu- 
ally cross one of the molecular energy level. Ac- 
cordingly, a current channel is opened and a jump 
in the I-V curve appears. The significant obser- 
vation is that, for the molecular bridge with BPD 
molecule (free from localizing group), the current 
amplitude is much higher (see solid curve of Fig. [S]) 
compared two the other two bridges. This is due to 
the fact that the localizing groups (both in BDMT 



and BPE molecules) interfere with the conjugated 
aromatic systems and suppress the overall conduc- 
tance through the molecules. On the other hand, 
the another important feature is that, in purely 
conjugate molecule (BPD) the electron conduction 
takes place as long as the bias voltage is applied, 
while for the other two molecules it appears beyond 
some finite values of V. This behavior gives a key 
idea in the fabrication of molecular devices. 

The shape and height of these current steps de- 
pend on the width of the molecular resonances. 
With the increase of molecule-to-electrode coupling 
strength, current gets a continuous variation with 
the applied bias voltage and achieves much higher 
values (compared to the current amplitude in the 
weak couphng case), as plotted in Fig.[6l where the 
solid, dotted and dashed curves correspond to the 
same meaning as in Fig. [5l In this strong molec- 
ular coupling limit, the current amplitude for the 
molecular bridge with BPE molecule is negligibly 
small compared to the other two bridges and the 
other features are also similar to the case of weak 
molecular coupling limit. 

4 Concluding remarks 

In summary, we have studied electron transport, at 
absolute zero temperature, through three short sin- 
gle conjugated molecules based on the tight-binding 
framework. We have used parametric approach, 
since we are interested only on the qualitative be- 
haviors instead of the quantitative ones, rather than 
the ab initio theories since the later theories are 
computationally too expensive. This technique can 
be used to study the electronic transport in any 
complicated molecular system. 

Electronic transport is significantly affected by 
(a) the molecule itself and (b) molecule-to-electrode 
coupling strength and in this article we have stud- 
ied our results in these aspects. In the weak- 
coupling limit conductance shows sharp resonant 
peaks, while these peaks get broadened in the limit 
of strong molecular coupling. These results predict 
that by tuning the molecular coupling strength one 
can control the electron conduction very sensitively 
through the molecular bridges. In the study of cur- 
rent we have seen that the current shows step-like 
behavior with sharp steps in the weak molecular 
coupling, while it becomes continuous in the strong- 
coupling limit as a function of applied bias voltage. 
Both for the two limiting cases our results have 
clearly described that the localizing groups suppress 
the current amplitude in large amount compared to 
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the current amplitude in case of purely conjugate 
molecule. Another significant observation is that 
the threshold bias voltage of electron conduction 
across a molecular bridge strongly depends on the 
molecule itself. These results provide key ideas for 
fabrication of different molecular devices, especially 
in the fabrication of molecular switches. 

Some assumptions have been taken into account 
for this present study. More studies are expected 
to take the Schottky effect which comes from the 
charge transfer across the molecule-electrode inter- 
faces, the static Stark effect, which is taken into 
account for the modification of the electronic struc- 
ture of the bridge system due to the applied bias 
voltage (essential especially for higher voltages). 
However, all these effects can be included into our 
framework by a simple generalization of the pre- 
sented formalism. In this article we have also ne- 
glected the effects of inelastic scattering processes 
and electron-electron correlation to characterize the 
electron transport through such bridges. 
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